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Abstract 

Voltage clamp data were analyzed in order to characterize the prop- 
erties of the fast potassium transient current Ia for a serotonergic neuron 
of the rat dorsal raphe nucleus (DRN) . We obtain maximal conductance, 
time constants of activation and inactivation, and the steady state acti- 
vation and inactivation functions m^ and hoc, as Boltzmann curves, de- 
fined by half-activation potentials and slope factors. We employ a novel 
jrt method to accurately obtain the activation function and compare the re- 

sults with those obtained by other methods. The form of Ia is estimated 
as g(V — V rev )m 4 h, with g — 20.5 nS. For activation, the half-activation 
potential is V a = —52.5 mV with slope factor k a = 16.5 mV, whereas for 
inactivation the corresponding quantities are -91.5 mV and -9.3 mV. We 
discuss the results in terms of the corresponding properties of Ia in other 
cell types and their possible relevance to pacemaking activity in 5-HT cells 
of the DRN. 

Keywords: Serotonin; dorsal raphe nucleus; potassium transient current; 
activation. 



1 Introduction 

Serotonergic neurons in the DRN, which extensively innervate most brain re- 
gions, have a large influence on many aspects of behavior, including sleep- wake 
cycles, mood and impulsivity (Liu ct al, 2004; Miyazaki et al., 2011). The fir- 
ing patterns of these cells have been much studied and many properties of the 
ionic currents underlying action potentials have been investigated (Aghajanian, 
1985; Segal, 1985; Burlhis and Aghajanian, 1987; Penington et al., 1991, 1992; 
Pcnington and Fox, 1995). In order to quantitatively analyze action potential 
generation in DRN 5-HT neurons, it is desirable to have accurate knowledge of 
the activation and inactivation properties of the various ion currents for some 
of which there is relatively little or no data available. Here we report results 
for the fast transient potassium current Ia which plays an important role in 
determining the cell's firing rate. 

Mathematical modeling of electrophysiological dynamics has been pursued 
for many different nerve and muscle cell types. Some well known neuronal 
examples are thalamic relay cells (Huguenard and McCormick, 1992; Destexhe 
et al., 1998; Rhodes and Llinas, 2005), dopaminergic cells (Komendantov ct 
al., 2004; Putzier et al., 2009; Kuznetsova et al., 2010), hippocampal pyramidal 
cells (Traub et al., 1991; Migliore et al., 1995; Poirazi et al., 2003; Xu and 
Clancy, 2008) and neocortical pyramidal cells (Destexhe et al., 2001; Traub 
ct al, 2003; Yu et al., 2008). Cardiac myocytes have also been the subject of 
numerous computational modeling studies with similar structure and equivalent 
complexity to that of neurons (Faber et al., 2007; Williams et al., 2010). 

The methods employed by Hodgkin and Huxley (1952) to describe mathe- 
matically the time (and space) course of nerve membrane potential have, for the 
most part, endured to the present day. An integral component of their model 
consists of differential equations for activation and (if present) inactivation vari- 
ables, gcncrically represented in the non-spatial models as m(t. V) and h(t, V), 
respectively, where t is time and V is membrane potential. These equations are 

dm mco — to dh ft-oo — h 
dt T m ' dt r h 

where TOoo(V^) and h^V) are steady state values and r m and t^ are time 
constants which often depend on V. For voltage-gated ion channels the current 
is often assumed to be 

I = g(V- V rev )m p h (2) 

where g is the maximal conductance, V rev is the reversal potential for the ion 
species under consideration and p is (usually) a small non-negative integer be- 
tween 1 and 4. The inactivation is invariably to the power one, as in (2). 
Although these equations have solutions which can only yield approximations 
to actual currents, it is nevertheless desirable to have accurate forms for the 
steady state activation and inactivation functions TOoo(T^) and /ioo(V), the time 
constants T m (V) and 77, (V), and the remaining parameters g and V rev . Our 
concern in this article is to describe and illustrate a novel and straighforward 
yet accurate method of estimating these quantities from voltage-clamp data. 



2 Results 

In order to isolate Ia, currents for activation and inactivation protocols were 
obtained in 20 niM TEA, as described in Section 4.1, for several identified sero- 
tonergic neurons of the rat dorsal raphe nucleus. In many cells the outward 
currents were clearly a mixture of components with different dynamical proper- 
ties. In some cells, however, the currents decayed smoothly from an early peak 
to near zero and were apparently a manifestation of a single channel type, pre- 
sumed to be uncontaminated Ia- Thus neither C0CI2 nor CdCl2 was employed. 
One cell in particular, DR5, with an apparently pure set of current traces was 
singled out for analysis. Its activation current traces are shown in Figure 1. Re- 
sults for other cells with composite outward currents will be analyzed in future 
articles. 




Figure 1: Activation current traces under voltage clamp for cell DR5 in 20 mM 
TEA. Time marker, 100 ms, Current marker, 100 pA. Voltage stepped from 
-120 mV to -60, -50, -40, -30 and -20 mV. For experimental protocol see Section 
4.1. 



The current traces of Figure 1 were obtained in digitized form and the time 
constants r m and Th and the value of p were estimated by the best fitting of the 
current traces to the analytical formula (6) . Only two of these three parameters 
can be independently chosen by virtue of the constraint given by formula (7) 



for the time of occurrence t max of the maximum current, which rearranges to 
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The current was required to pass through the maximum point at (t max ,I max ) 
and one other point on the trace, called (£2, Ii)- The value of p = 4 was found 
to give excellent fits to the current traces. An example of the resulting curve 
with best fitting time constants r m and Tu is shown in Figure 2 along with the 
experimental current for a step from -120 mV to -20 mV. 

Table 1 summarizes the data and estimates of time constants for cell DR5, 
all of which quantities are required to estimate the three parameters g (conduc- 
tance), V a and k a (Boltzmann parameters for activation). In this table Ia, max 
is the maximal experimental current over time. 
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Figure 2: An example of the best fitting curve for the calculated time-dependent 
current Ia(£) (smooth curve), as given in Equation (6) with p = 4, to an experi- 
mental trace for cell DR5 with calculated time constants for activation r m = 1.5 
ms, and inactivation r^ = 28 ms. Voltage stepped from -120 mV to -20 mV. 



The results obtained by nonlinear least-squares fitting (Section 4.2.2) with 
correction factor as in Equation (8) and following but without making any 
assumption about V* , (called method A), are that for cell DR5 the whole cell 
maximal conductance for the Ia channels is g — 20.5 nS, the half-activation 



Table 1: Maximal currents and estimated time constants in 
activation experiments for Ia in cell DR5 

Step from -120 mV lA,max (pA) r m (ms) r h (ms) 
to 



-20 mV 


825.4 


1.5 


28.0 


-30 mV 


431.7 


1.5 


28.0 


-40 mV 


171.5 


2.4 


21.7 


-50 mV 


22.2 


- 


- 


-60 mV 





- 


- 



potential is V a — —52.5 mV and the corresponding slope factor is k a = 16.5 mV 
(see also Table 2). 

As a check on the estimated parameter values by method A, the maximal 
currents with steps from Vq = —120 mV to V\ =-60, -50, -40, -30 and -20 mV 
were calculated by formula (8) using the values g = 20.5 nS, V a = —52.5 mV, 
k a = 16.5 mV, p = 4, V rev = —105 mV and the time constants given in Table 1. 
The maximal currents so determined arc plotted against V\ in Figure 3 where 
it can be seen that there is excellent agreement with the experimental values. 
Also shown in the inset of the Figure is the activation function raised to the 
power 4, a so called 4th order Boltzmann. 
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It is interesting to compare the values of the activation function so obtained 



Figure 3: Predicted and experimental values of I m ax in activation protocols 
stepped up from -120 mV to various clamp voltages, in mV, -60, -50, -40, -30 
and -20, for cell DR5 with TEA 20 mM. The inset shows the best fitting 4th 
order Boltzmann. 

with those obtained by three other procedures. For method B, in which it is 
assumed that V* = —20 mV so that moo (—20) = 1 (see note after Equation 
(11)), which implies that Equation (12) gives an estimate of g, the values g = 
12.91 nS, V a = -54.7 mV and k a = 12.52 mV are obtained. 

In method C, the value of V* is not fixed but a best-fitting Boltzmann is 
obtained by least sqares, which gives V a = —47.0 mV and k a = 10.2 mV. It is 
also of interest to compare these sets of results with those by the usual method, 
called method D, using the raw G/G max values which assumes that V* = — 20 
mV. This gives V a = -48.2 mV and k a = 11.78 mV. 

It can be seen from the data in Table 2 that there are considerable discrep- 
ancies between the results obtained by the various methods. Method A gives 
the most accurate estimates. For the other three methods, relative errors in 
the conductance are between -37.1 % and -52.6 % relative errors in V a are from 
-10.5 % to +4.2 % and those in k a are from -24.1 % to -38.2 %. Such errors can 
clearly lead to significant errors in the computed ionic currents. 



Table 2: Estimated activation function and conductance of Ia 

for cell DR5, obtained by four methods as described in the 

text: % errors relative to method A given in brackets 



Method -> 


A 


B 


C 


D 


5(nS) 


20.5 


12.9 (-37.1) 


12.4 (-39.5) 


9.71 (-52.6) 


V a (mV) 


-52.5 


-54.7 (+4.2) 


-47 (-10.5) 


-48.2 (-8.2) 


k a (mV) 


16.5 


12.52 (-24.1) 


10.2 (-38.2) 


11.78 (-28.6) 



2.1 Inactivation experiments 

For the inactivation voltage-clamp experiments, whose results are not shown, the 
maximum currents J maa; (Vo, V\) at various starting potentials Vo are sufficient 
to estimate the relative values of the inactivation function because all factors 
on the right hand side of (17) except /loo(Vo) do not depend on Vq. Assuming 
that hoc(— 120) = 1, the best fitting steady state inactivation function was 

h °°(V) = 1 + e (V + 91.5)/9.3 - W 

2.2 Graphical summary 

The activation functions obtained by methods A-D and the inactivation function 
for the transient potassium current Ia in cell DR5 are drawn in the top panel of 



Figure 4. Also shown, in the lower panel, is the activation function to power 4, 
which reflects more closely the contribution of activation to current amplitudes 
at various membrane potentials. 
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Figure 4: Top part. Steady state inactivation function and steady state activa- 
tion function obtained by the methods A-D as explained in the text for Ia in cell 
DR5 obtained from voltage-clamp data by the methods described. The result 
for method A is the most accurate estimate of ra^ . Lower part. The activation 
function to the power 4 is shown for the methods A-D. 



3 Discussion 



There have been many investigations of the properties of the ionic currents 
involved in the firing of serotonergic neurons of the dorsal raphe nucleus (Agha- 
janian, 1985; Segal, 1985; Burlhis and Aghajanian, 1987; Penington et al., 1991, 
1992; Penington and Fox, 1995; Chen et al., 2002). However, the basis of appar- 
ent pacemaker-like activity in these cells is not fully understood. It was claimed 
that noradrenergic input is required (Burlhis and Aghajanian, 1987) for firing 
in some cells and there have been conflicting ideas about the role of Ia (Burlhis 
and Aghajanian, 1987; Segal, 1985). No reports on the properties of Ia in these 
cells have previously appeared. 



In the first quantitative study of Ia, Connor and Stevens (1971) reported 
an m A h form for the current in a single gastropod neuron and found quite 
long time constants for activation (12 ms) and inactivation (235 ms). The V1/2 
value for inactivation was -68 mV, with a resting potential of -40 mV. It was 
also found that Ia was important mainly in the middle and end of the inter- 
spike interval (ISI, see below). Concerning the method of analyzing voltage- 
clamp data for activation, in their pioneering work on squid axon, for large 
depolarizations Hodgkin and Huxley (1952) used the general Equation (5) for 
the time-dependent sodium current converted to a conductance 

g(t) = ff[moo(Vi)(l - e- t / T -)] 3 / loo (Fo)e- t/T " 

to estimate 5[m o(Vi)] 3 ft- 00 (Vo). Bclluzzi and Sacchi (1988) also employed Equa- 
tion (5) with simplifying assumptions for m and h, and pas a parameter. Ac- 
tivation time constants were estimated via Equation (7) with p — 3, the small 
values of both r m and Th they obtained reflecting the relatively high temperature 
of the preparation. For activation in histamine neurons, Greene et al. (1990) 
found a fast and slow component of Ia and reported Boltzmann parameters for 
normalized current without converting to conductances. Huguenard and Mc- 
Cormick (1992) used normalized currents to find inactivation parameters and 
time-courses of conductances (as in (5)) to determine Boltzmann parameters for 
activation. Bekkers (2000) gave two sets of results for activation, one obtained 
from the empirical Boltzmann fit to the peak Ia conductance curve, and the 
other obtained from the time-dependent Ia conductance curves. We have pur- 
sued a new approach, as outlined in the last section, which utilizes the maxima 
of the currents at various clamp voltages and employed a correction factor as in 
Equation (8). 

The results obtained for Ia in a serotonergic cell of rat DRN are compared 
with those in some other CNS preparations in Table 3. It can be seen that 
the half-activation potentials for activation are generally similar and that the 
largest variation occurs in the time constant of inactivation. 

In order to construct a computational model of action potential generation 
in DRN 5-HT neurons, such that the role of the various ionic components can be 
properly understood, it is desirable to have accurate knowledge of the activation 
and inactivation properties of the various ion currents, including that addressed 
here, Ia- The need for considerable accuracy is made the more important be- 
cause several components, including the low threshold T-type calcium and Ia 
operate in overlapping ranges of potential. In this paper we have illustrated 
that the often-used procedure of finding activation functions ntoo(Vi) by taking 
ratios of current to maximal current or conductance to maximal conductance 
may produce inaccurate results due to the fact that the quantity F p {^f) in (8) 
depends on the clamp voltage and does not cancel. Not taking this into account 
can lead to substantial errors in the estimates of the parameters k a and V a . 
Inactivation functions, however, do not have such a complication and can be 
estimated in the usual way. 



Table 3: Activation and inactivation parameters for Ia in some CNS cells 





DRN 

serotonergic" 


Sympathetic 
cervical ganglion 6 


Thalamic 
relay c 


Cortical 
pyramid^ 


Activation 


V1/2 mV 
k mV 
T m ms 


-52.5 

16.5 

1.5-2.4 


-58.55 

14.39 

< 1 


-60 

8.5 

0.5-2.5 


-56.2 
19.1 
0.3-8 



Inactivation 










Vi/2 mV 
k mV 

T h ms 


-91.5 
9.3 

21.7-28.0 


-78 

7.3 

< 10 


-78 

6 

12-65 


-81.6 
6.7 
6-8 



"This study; ^Belluzzi & Sacchi (1988); c Huguenard & McCormick (1992); 

d Bekkers (2000). 



L-type calcium currents, particularly through Ca v 1.3 channels, have been 
found to sometimes play a role in pacemaker activity in neurons and cardiac cells 
(Koschak et al., 2003; Putzier et al., 2009; Marcantoni et al., 2010). However, in 
DRN 5-HT cells, there is only a small (about 4% of total lea) contribution from 
L-type calcium currents (Penington et al., 1991) so it is yet to be determined 
what role these play in the pacemaking activity. These latter results were, how- 
ever, obtained in dissociated cells. Burlhis and Aghajanian (1987) hypothesized 
that the low threshold calcium T-type current played a key role in pacemaking 
and Segal (1985) proposed that Ia was also important. However, considering 
the fact that the resting potential for these cells is about -60 mV and that 
threshold for spiking is about -50 mV, and judging by the curves for to^V) 
and hoo(V) in Figure 4, it is likely that Ia is only important during the spike 
itself and not between the commencement of the afterhyperpolarization and the 
next spike, probably implying a lesser role for Ia in pacemaking as it tends to 
be switched off during the greater part of the ISI. We will explore by means of a 
computational model the ionic basis of the mechanisms of pacemaking in DRN 
5-HT cells in future articles. 

Finally we note that the Ia results obtained here for m^ and the time 
constants bear a resemblance to those given in Gutman et al. (2005) for the 
transient A-type channels K v 4.1 ( V a = —47.9 mV, k a — 24 mV, Th — 22 ms, 
and two other Th values). This contrasts with the finding by Bckkers (2000) 
who suggested that in rat layer 5 cortical pyramidal cells, Ia was carried by 
K V A.2 channels, which accounts for their smaller time constants of inactivation 
(Gutman et al., 2005). 



4 Materials and Methods 

4.1 Experimental 

Data of the kind analyzed in this paper were obtained from 2 month-old male 
Spraguc-Dawley rats that were anesthetized with halothane and then decapi- 
tated with a small animal guillotine in accordance with our local Animal Care 
and Use Committee regulations. Three coronal slices (500 /zm) through the 
brain stem at the level of the DRN were prepared using a vibroslicc in a manner 
that has previously been described (Yao et al., 2010). Pieces of tissue contain- 
ing the DRN were then incubated in a PIPES buffer solution containing 0.2 
mg/mL trypsin (Sigma Type XI) under pure oxygen for 120 min. The tissue 
was then triturated in Dulbecco's modified Eagle's medium to yield individual 
neuronal cell bodies with vestigial processes. Recording was carried out at room 
temperature, 20°C-25°C. Cells recorded from had diameters of about 20 /jm. 
Stcinbusch ct al. (1981) showed using immunohistochcmistry that most of the 
cells in a thin raphe slice with a soma diameter greater than 20 /im contain 5-HT 
while the smaller cells were largely GABAergic interneurons. Using a method 
similar to that of Yamamoto et al. (1981), we also found that the proportion of 
isolated cells with diameters larger than 20 /j,m that stain for 5-HT, was greater 
than 85%, and a similar percentage responded to serotonin. 

The cells were voltage-clamped using a switching voltage clamp amplifier 
(Axoclamp 2A) and a single patch pipette in the whole-cell configuration (Hamill 
et al., 1981). Pipettes pulled from thick-walled borosilicate glass (resistance 
from 5-7 MSI) allowed a switching frequency of 8-13 KHz with a 30% duty 
cycle. The seal resistance measured by the voltage response to a 50 pA step of 
current was often greater than 5 Gil). The voltage-clamp data were filtered at 
1 KHz and digitized for storage at 16 bits, or experiments were run online with 
a PC and a CED 1401 interface. 

The external saline was designed to isolate potassium currents and contained: 
outside the cell TTX 0.2 (M and in mM NaCl 147, KC1 2.5, CaCl 2 2, MgCl 2 2, 
Glucose 10, HEPES 20, ph 7.3. Inside: K+ Gluconate 84, MgATP 2, KC1 38, 
EGTA-KOH 11, HEPES 10, CaCl 2 1, pH 7.3. The total [K+] inside was 155mM 
(since 33 mM is added to make the salt of EGTA). The osmolarity was adjusted 
with sucrose so that the pipette solution was 10 mOsm hypoosmotic to the bath 
solution. Drugs were either dissolved in the extracellular solution and added to 
the perfusate or applied by diffusion from a patch pipette (tip 15 /im) lowered 
close to the cell (50 /im and then removed from the bath). As a control, the 
same procedure was carried out without addition of the drug to the application 
pipette; this did not affect Ia- In some cells, to investigate any influence of 
Ca 2+ currents on Ia either 2 mM C0CI2 or 20 /jM CdCl 2 was added to the bath 
by a pipette placed near the cell. This sometimes changed the magnitude of I a 
but had little effect on its time course. 
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4.2 Theory 

Suppose that in a voltage clamp experiment the voltage is, after equilibration 
at a voltage Vo, suddenly clamped at the new voltage V\. Then, assuming a 
current of the form of Equation (2), according to the standard Hodgkin-Huxlcy 
(1952) theory, the current at time t after the switch to V\ is 

/(*; Vo, Vi) - g(V - V rev )[ mi - (mi - m Q )e- t/T ™] p [hi - (hi - Me - ' 7 "*], (5) 

where we have employed the abbreviations too = m oc (Vo) 1 mi — TOoo(Vi) for 
the activation, ho = /i<x>(Vo),/ii = hoo{Vi) for the inactivation and r m = 
T m (Vi), Th — Th{Vi) for the time constants. 

4.2.1 Activation experiments 

In activation eperiments there is a step up from a relatively hyperpolarized 
state Vo to several more depolarized states V\ so that one usually assumes that 
ho — 1, hi — and too = 0. This gives the simplified expression for the current 

I(t; Vi) = g(V - V rev )[ moo (Vi)(l - e^™)]^-^* . (6) 

Finding the maximum by differentiation yields the time at which the maximum 
occurs as 

t m ax{Vl)=T m \u(l + ^] (7) 



and its value as 
where 



Imax{Vl) = g(Vi - Vr^mUV^F^) (8) 

F ^ - (TT^ (9) 



and where 7, which depends on V\, is defined as 

7(*) = ^£l (10) 

T m ( Vl ) 

Since V\ is here considered a variable, the time course of the current changes as 
Vi varies. By rearrangement the steady state activation is given by 

tir \ I J-max(Vl) \ p fti\ 

m - {Vl) = Kwi-v^m,)) (11) 

If it is known that for some (usually well-depolarized) state V — V* one has 
Trioo{V*) = 1 (but note that this can only ever be approximately true as TOqo 
can only approach unity asymptotically) then the maximal conductance g can 
be estimated from voltage clamp experiments from 

77= imaxiV*) ( 2] 

9 (V*-V rev )F p ( 7 (V*)Y { ] 
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Once g is known then (11) can be used to find moo(Vi) for various V\ because 
all remaining quantities, I maxi (V\ — V rev ) and F p (-y(Vi)) are known. 
In summary, the steps to find the activation function moo(l^) are: 

(a) From the time course of the current find the maximum current and its time 
of occurrence for a given V\ 

(b) Estimate the time constants T m {V\) and r^(Vi) and the value of p. 

(c) Hence find 7(Vi) and also F p (j(Vi)). 

(d) Using a suitably highly depolarized state V — V* estimate g according to 

(12). 

(e) Hence estimate m^^) at various V below V* . 

(f ) Find the best fitting Boltzmann curve passing through the values obtained. 
However, the formulas (6)-(8) were derived on the assumption that h{V\) = 1 

and this will usually be valid for V c < V\ < V* , where V c is below the half- 
activation potential V a . This means that the half-activation function moo(V) 
can be obtained accurately for V > V a , and thus by symmetry for V < V a . 
Another approach is as follows. 

4.2.2 V* is not known: least squares estimates 

If V* is not known or is uncertain, then least squares estimates can be made as 
follows. It is assumed that 

m -( y ) = i + e -(v-v a) /k a ( 13 ) 

is the steady state activation function. Then one may estimate g, V a and k a 
by minimizing the sum of squares of the differences between experimentally 
obtained maxima of the current at various Vi and the values given by formula 

(8) 

4.2.3 Summary of methods A-D used to estimate the activation 
function 

Here is given a brief description of the 4 methods employed to estimate the 
parameters of m^V) and g given in Table 2. 

Method A 

Having estimated the values of the time constants r TO (Vj) and Th(Vi) and the 
power p as described in the text, experimental maximum currents at various Vi 
are compared with those obtained from formula (8) 

I max (Vi) - g(Vi - Vr^m^WF^) (14) 

with 7 evaluated at Vi. Thus the parameters g, V a and k a are estimated by 
least squares. 

Method B 

It is assumed that moo(V*) — 1, with V* = —20 mV in the present example, 
so that g is given directly by Equation (12). Then m^Vi) can be obtained for 
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various Vi uising Equation (11). A best fitting Boltzmann is then obtained for 
the various vaules of m^Vi). 

Method C 

The experimental values of the maximum currents (over time) Ii. m ax(Vi) at 
each test voltage V, are divided by Vi — V rev to give corresponding conductances 
Gi(Vi). The conductances Gi(V) are assumed (from Equation (8) without the 
correction factor) to be given by 

G i (V i )=gm* (V i ). (15) 

With 

m -( y ) = i + e -(v-v a) /k a ( 16 ) 

the parameters ~g, V a and k a (and possibly p if not already known) are estimated 
by least squares fitting of the experimental Gi values to the predicted ones, 
without any assumption on V* (normalization). 

Method D 

The same procedure is carried out as in method C, but the Gi values are divided 
by the maximum value G max , assumed to occur at V* (-20 mV in the above), 
to give the normalized values 

G, = -^-. (17) 

^max 

Then since Gi (V* ) = 1 there are only two parameters V a and k a to estimate with 
a best fitting Boltzmann curve. The maximal conductance g can be obtained 
as 

* * rev 

4.2.4 No inactivation 

If there is no inactivation in the course of the clamp experiment then the max- 
iumum current will be the asymptotic value 

I max (Vi) = g(Vi ~ ^XW). (19) 

So, with V* as above, g is again given by (18), and it is straightforward to 
obtain the required activation function from the remaining I max measurements. 

4.2.5 Inactivation experiments 

Here steps are made from a number of relatively hyperpolarized states Vq to 
a fixed relatively depolarized state V\. It is usually assumed that h\ — and 
rriQ = 0, so that the current is approximately 

I(t; V , V x ) = g(Vt - F re ,)mS (F 1 )/i (X) (Fo)(l - e^™)^-*/ 7 *. (20) 
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Here the time course of the current (but not its magnitude) is the same for all 
V . 

The time of occurrence of the maximum t max of I(t; Vq; V\) is again given 
by (7) but the value of the maximum is now 

ImaAVo, V x ) = g{Vi - Vr^mUVi^^Fp^). (21) 

Assuming that g has already been estimated and that V* is such that 
m-ooiV*) = 1, then the inactivation function can be found at Vq from 

hoo(V ) = I max (Vo, V*)g(V* - U re „)F p ( 7 (U*)) (22) 

It can be seen that normalizing the currents by dividing by the maximum of the 
I max gives a reasonable estimate of the steady state inactivation function. 
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